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Abstract 

Recent experimental and theoretical progress regarding B — > K^l + l~ de- 
cays led to improved bounds on the Wilson coefficients Cg and C\q of four- 
fermion operators of the |A£?| = |A«S| = 1 effective Hamiltonian. We an- 
alyze the resulting implications on squark flavor violation in the MSSM and 
obtain new constraints on flavor-changing left-right mixing in the up-squark- 
sector. We find the dimensionless flavor mixing parameter (S^lr, depending 
on the flavor-diagonal MSSM masses and couplings, to be as low as < 0.1. 
This has implications for models based on radiative flavor violation and leads 
to B{B S —7- > 1 x 1CP 9 . Rare top decays t — > cj,t — > cg,t — > cZ have 

branching ratios predicted to be below < few x 10~ 8 , 1CT 6 and 10 -7 , respectively. 



1 Introduction 



Heavy flavor physics is rapidly advancing with the successful start of the Large Hadron 
Collider's (LHC) 6-physics program and the final analyses from the Tevatron as well 
as the i?-factory experiments Belle and BaBar. Most notably the LHCb collaboration 
is currently making a clean sweep in model space around the Standard Model (SM): 
The current upper limit on the B s fi+fT branching ratio at 95% (90%) C.L. (I) 

B(B S ->■ < 4.5 (3.8) x HT 9 (1) 



is down to the level of the SM, B(B S / u + / u") SM ~ (2.5 - 4) x 10 j 3 , e.g., |2||3j. The 
SM prediction is much more precise assuming the measured B s — B s mass difference 
to be SM-like, B(B S -> // + //~)sm = (3.2 ± 0.2) x 10~ 9 |4|. Another milestone consti- 
tutes the preliminary measurement of the position of the zero of the forward-backward 
asymmetry in B — > K°* // + decays |5| 



g 2 =4.9^GeV 2 , (2) 

consistent with the SM prediction g 2 | SM = 4.0 ± 0.3 GeV 2 §, also [7|[8). 

In this work we aim at investigating the space for supersymmetric flavor physics 
in the light of the recent and new data from direct collider searches and on flavor- 
changing rare processes b — > Especially, the study of the exclusive modes 
B — > K^fi + fi~, which are accessible to hadron colliders, has progressed significantly 
over the last year(s) both experimentally and theoretically. The latter is due to the 
exploitation of the region where the invariant mass of the dilepton is large, of the 
order of the 6-quark mass |9|, by making use of the heavy quark effective theory 



framework of Ref. 10 . Using the recent data from CDF 11 and LHCb 12 , improved 
model-independent constraints on the Wilson coefficients Cg^o of the semileptonic four- 
fermion operators 9 and O w have been obtained (6,13 



Here, we study the implications of these Cg^o-constraints on squark flavor viola- 
tion in the minimal supersymmetric SM (MSSM). Flavor violation in supersymmetric 
(SUSY) models originates from generational mixings in the sfermion mass matrices. 
One-loop SUSY effects including their phenomenology in semileptonic 5-decays are 



known for a while 14-17 . We work out the implications of the improved data on 
semileptonic decays, taking into account the recent direct search limits on SUSY par- 
ticles. 

The plan of this paper is as follows: In Section [2] we explain, after introducing 
the | AS | = 1 effective theory framework, why only chirality-flipping flavor mixing 
between the second and third generation in the up-sector, parametrized as (S^lr, has 
a significant sensitivity to the C 9> io-constramts. We then describe the method of SUSY 
parameter scanning and discuss the constraints on (5% 3 ) lr- Results are given in Section 
|3j including predictions for &-physics, implications for the models with radiative flavor 



violation (RFV) recently discussed in Refs. 18 19 , and rare top decays. We summarize 
in Section |4| Details on the SUSY loop contributions to b — > sl + l~ processes including 
corrections to the literature are given in the Appendix. 
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2 Constraining squark flavor mixing 



In Section 2.1 the effective theory for \AB\ = 1 transitions is described. In Section 2.2 



we discuss which flavor parameters in supersymmetric models can be effectively probed 



in semileptonic decays. The constraints are worked out in Section |2.3[ and discussed 
in Section EH 

2.1 The \AB\ = 1 effective theory framework 

In order to describe rare decays of 5-mesons, we employ the effective Hamiltonian 

AG, 



H 



cff 



~^V tb V* CiMOifr) + h.c. , (3) 



where the Cj are Wilson coefficients and the Oj higher-dimensional \AB\ = \AS\ = 1 
operators, and V denotes the Cabibbo-Kobayashi-Maskawa (CKM) matrix, Gp the 
Fermi constant and fi the factorization scale. The New Physics (NP) contribution 
to the current-current (Oi^), QCD penguin (O^-q) and chromomagnetic dipole (Og) 
operators can be neglected here (cf., e.g., |14|). The most relevant operators are the 
electromagnetic dipole operator 7 and the semileptonic four-fermion operators 9tW , 
written as 

9 = (s L l,b L ) (m) , Oio = ^ {HlM (h^l) ■ (4) 

The MS mass of the 6-quark is denoted by m&, and we neglect the mass of the strange 
quark. 

In general models for NP, the effective Hamiltonian contains additional operators 
with flipped chirality. In the SM and in minimally flavor-violating (MFV) SUSY mod- 
els with CKM-induced flavor violation, these are suppressed by m s /rrib compared to 
the corresponding unflipped operators given in Eq. Q. The non-MFV chargino con- 
tributions to b sl + l~ processes share this feature of suppressed chirality-flipped con- 
tributions. As discussed below, we will be mainly interested in chargino loop-induced 
contributions, so it is justified to neglect the chirality-flipped operators and consider 
left-handed currents only. Furthermore, we do not take into account the effects of the 
scalar- and pseudoscalar semileptonic operators (and their chirality-flipped counter- 



parts). This is a good approximation for not too large tan/3 20 , to which we restrict 
our analysis. The exclusion of a sizable enhancement in the B s — > branching 
ratio, see Eq. (JTJ, supports this further. 

The NP contribution, which in our case is the MSSM contribution, can be split into 
one stemming from the diagonal elements of the squark mass matrices, labeled 'diag', 
and the remainder induced by the corresponding flavor off-diagonal entries as 

Cf p = Cf ag + 1 , (5) 
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173.3 GeV 


25 




m h {m b ) 


4.19 GeV 


26 




m w 


80.399 GeV 


26 




m z 


91.1876 GeV 


26 




a s {m z ) 
s 2 


0.1184 


26 




0.23116 


26 





Table 1: The values of the Standard Model parameters used in this work. 



where Cf p denotes the NP contribution to the C u i.e., d = Cf M + Cf p . We use the 
terminology 'mass insertion' (MI) since we will constrain the commonly used MI pa- 
rameters, i.e., the off-diagonal elements of the squark mass term divided by an average 
squark mass squared. Note, however, that we do not rely on the MI approximation 



21 , 22 in the numerical analysis, but instead use the exact formulae. 



We employ the two-loop matching conditions of Ref. 23 for the SM contribution 



and the one-loop results of 14 16 for the MSSM contribution. The Wilson coefficients 
at the relevant scale for 6-decays, = 0(mb), are obtained from the ones at the 
matching scale /i of (9(100 GeV) by solving the RG equations. We perform this by 
extending and using the flavor tool EOS 



24 



The values of the SM parameters used 
are given in Table [TJ The SM values of the most important Wilson coefficients are 
given as (at ^ = 4.2 GeV) 



CfV) = -0.33, 



CfV & ) = 4.27, 



Cfo M M = -4-15. 



(6) 



2.2 SUSY flavor contributions to b — > sl + l 

We begin by discussing which MI parameters have the best potential of receiving 
significantly improved constraints from the Cg^o-bounds. Relevant for Cg^o are (S^ll 
and (#23) lr, which enter through gluino loops, and (5% 3 ) ll and (^3) lr, which appear 
in chargino loops. Of these, the ones which are not yet very much constrained by 
bounds on other Wilson coefficients (C7 in this case) and at the same time give a 
substantial contribution to Cg,io are most interesting: 

• The parameter (5$ 3 ) lr is already tightly constrained by bounds on C 7 from data 
on the B — > X s ^ branching ratio 22 . In addition, (S^lr contributes to Cg,io 



only in double MI diagrams [16). We conclude that it plays no role for our 
analysis. 

The parameter (5| 3 )ll is much less constrained by B — > X s ^ and B s — B s mixing 



27 



than (^23) lr. On the other hand, (^23)^ has little effect on Cg^o: Its 
effect in the Z-penguin gluino loop is suppressed with respect to that in the 
7-penguin by a factor of m 2 /m| (as discussed, e.g., in [16]). The 7-penguin in 
turn does not contribute at all to C w , and its contribution to Cg is numerically 
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tan/3 


M 2 






rrig 






m 9 


300 


4 


150 


-300 


300 


1000 


1000 


100 


700 



Table 2: Example SUSY point at = 120 GeV. All masses are in GeV. m§, are 
the diagonal elements of the squark mass matrix as given in the Appendix. The point 
corresponds to the following spectrum: m% = 236 GeV, m ( - 2 = 1017 GeV, = 150 
GeV, mc 2 = 321 GeV and m h o = 117 GeV. 





qT^o) 




0M(5% 3 ) LR -0.38(5% 3 ) LL 


0.17(^3)^-0.11(^3)^ 


-2.24(5« 3 ) Lfi + 0.19(5 2 « 3 ) LL 




GT^o) 


C(w) 


16.35 {5 d 23 ) LR -0.02(^ 3 ) ll 


0.04 ($ 3 ) LL 





Table 3: Flavor off-diagonal Wilson coefficients at /to = 120 GeV for the SUSY point 
given in Table [2j Double mass insertions are not shown. 



subleading to those from chargino loops, in particular for squark masses above a 
TeV. Therefore, the parameter (S^ll is not of interest to our study. 

• The coefficients (#23)1,1, an d (S^lr both give a non-negligible contribution to 
C7, which is, however, by far not as significant as that from (S$ 3 )lr- The impact 
of ($33) LX-, X = L,R contributions to B s — B s mixing is very small and negligible 
[27| . As exemplified in Table [3] for the SUSY benchmark point given in Table § 
(&2z)ll gives a much larger contribution to C 7 than (S^lr- Moreover, C% is 
about an order of magnitude more sensitive to (8% 3 )lr than to (S^ll- We 
conclude that (#23) lr, is the most relevant parameter for our analysis on updated 
constraints from C^io-bounds. 

For completeness and because of disagreement with the literature in the photon- 
penguin [l6j, see Appendix, we give here the dependence of (77,9,10 on (5% 3 )lr in the 
MI approximation, 

yMi,x{ \ _ K*s r^w 



Cr' X (»o) = ^-^FxW 3 ) LR 
V ts 92 mq 



n Ml ^(n \ cs 1 / ,2 1 \ TpZ-p , .2 171 W F 7-p %pboxV AM 

C 9 (MoJ - ~U 4s W - + %^2" )l d 23)M 



'*S ^^VV ^2 - ""g 



where At denotes the top Yukawa coupling and the sine squared of the weak mixing 
angle. The expressions F and _F 7_P are due to 7-penguin-diagrams, while F z ~ p and 
F box are due to Z-penguin- and box-diagrams, respectively. Their explicit form can be 
found in the Appendix. The up-squarks are assumed to be roughly degenerate, except 
for the right-handed stop, which is allowed to be significantly lighter. Notice that 
the 7-penguin and box-contributions are parametrically suppressed by m^/m| with 
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respect to the Z-penguin effect. Furthermore, since the Z-boson mainly has an axial- 
vector coupling to charged leptons, the Z-penguin contribution to C 9 is suppressed to 
the one to C w by \As 2 w - 1| < 1 (28). 



2.3 Numerical analysis 

For the numerical analysis we scan the SUSY parameter space within the ranges given 
in Table |4} All parameters are assumed to be at the electroweak scale, taken to 
be Hq = 120 GeV. We fix the lightest sneutrino mass rriy = 100 GeV because the 
dependence of C^f on mc is very mild only. We denote by m~ the kth, k — 1, . . . , 5 
diagonal element of the up-squark mass matrix and by m~ the remaining diagonal 
entry connected to the right-handed stop, including the F and D terms, see Appendix. 
We further fix rrig = 1000 GeV because the bounds can get only weaker for larger 
average squark masses. All flavor- violating MI parameters except for (#2 3 )lr are set 
to zero. 

Then, for each parameter point we require the following constraints to be satisfied: 

i) The bounds from the B — > X s j branching ratio: 0.3 < \C 7 (fib)\ < 0.4 [6j, see 
text. 



it) The lightest chargino mass limit: m x ± > 94 GeV [26[ . 



Hi) The lightest stop mass limit > 100 GeV from D0 29 . We choose this lower 
bound since the stronger one from CDF, > 180 GeV [30], assumes the stop 
predominantly decaying via charged currents to 6-quarks, a lepton and missing 



energy. The recent LHC findings 31,32 are model-dependent, too, and not taken 
into account here. 



iv) The Higgs mass limits: m^o > 114.4 GeV 26 and m^o > 200 GeV. We require 
the latter to ensure that h° is sufficiently SM-like, see text. 



v) Electroweak precision tests: -0.0007 < Ap < 0.0017 26 . 





tan f3 


m H ± 


M 2 


H 


m ~t H 


A t 




min. 


3 


300 


100 


80 


170 


-3000 


-0.85 


max. 


15 


1000 


1000 


1000 


800 


3000 


0.85 



Table 4: Ranges of the SUSY parameters (at /io = 120 GeV) used for the parameter 
scan, see text. All masses are in GeV. We fix = 100 GeV and rrig = 1000 GeV. 



The MSSM allows for large contributions to CV, which can flip its sign while still 
being in agreement with data, see i). However, the existence of a zero-crossing of the 
forward-backward asymmetry in B — > K*fi + fi~, see Eq. enforces sign(CVCg) to be 
SM-like. Since C$ cannot change sign within the MSSM |17| , we allow in our analysis 
only MSSM points with SM-like signed C 7 . 
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We use the FeynHiggs code 33-36 v.2.9.0-betsQto check iv) and v). Further, we 
add 3 GeV to the prediction for m^o in order to account for the theoretical uncertainty 
[35] . Note further that the recent 95% CL. exclusion limits for a SM-like Higgs boson 
from ATLAS and CMS suggest a valid range around (125 ± few) GeV [37|B- We 
checked that while such values require typically larger A t than the constraint iv), the 
generic features of our flavor analysis hold. 

— 2 

Throughout this work we choose the normalization (5% 3 )lr = (A% 3 ) LR /M U , where 
M u = (f^rc.- + \ m 'i R ) to connect the dimensionless MI parameter (5% 3 )lr to the off- 
diagonal element {^23) lr of the up-squark mass matrix in the super-CKM basis. Note 
that values of (£33) lr outside the range given in Table H would yield tachyonic squarks 
due to the large, order one off-diagonal entry. Note also that a normalization to the 
geometric mean of the diagonal entries would result in increasingly larger values of 
(#23) lr for lighter right-handed stops mf < Trig. 

We begin by assuming flavor-diagonal SUSY contributions only, i.e., for (5% 3 )lr = 0. 
In this case charged Higgs boson and CKM-driven chargino loops give non-vanishing 
NP contributions to the Wilson coefficients. Defining the ratios 



Ri 



a 



NP 



c SM 

we find the following, maximally possible NP effects 

R*M<3%, RM<11%. (MFV) (9) 



As expected [17], the flavor-diagonal SUSY effects on Cg^o are small and well within 
what is allowed by present and near-future data. 

Switching on squark flavor violation, the possible size of NP effects in the Wilson 
coefficients increases 

RM<4:%, R w {n b )<47%, (no s.l. bounds) (10) 

in particular in C\q. Here, the new semileptonic bounds are not yet taken into account. 



Note that the Higgs mass bound iv) is efficient, as previously noted in 39 ; if ignored, 
the allowed range for Rio(fib) would be more than a factor of two larger. 

The outcome of the full scan is shown in Fig. [T] SUSY points that pass i) - v) are 
shown together with the allowed regions (gray areas) in the Cg(/ib) — Cio(^b) plane 
from the recent analysis of B — > K^*H + l~ decays [6] for SM-like signed C7. The SUSY 
flavor effects are dominated by the Z-penguin contribution, correlating Cg and Ciq as 
Cio~ P I Cg~ p = l/(4spy — 1), see Eq. ([7j). The latter correlation is shown as the dotted 
line and is clearly featured by the model. The analogous solutions with flipped-sign 
C-j > 0, which are allowed model-independently |6), are excluded in the MSSM (and 
not shown): firstly because the Z-penguin dominated scatter points miss in this case 
the allowed Cg(/x&) — Cio(//b) parameter space and secondly because of the measurement 
Eq. (g. 

: We are grateful to Sven Heinemeyer for FeynHiggs support and provision of version 2.9.0-beta. 
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Figure 1: Reach of SUSY models with (S^lr 7^ in the CgO^^Cio (/•*&) plane for 
C 7 (fj,b) < 0. The light (dark) gray shaded areas are the 95% (68%) confidence limit 
bounds which were obtained from B — > K^*H + l~ data in Fig. 7 of Ref. jfj]. The red 
dotted line denotes the Z-penguin correlation Cf ~ P /Cg p = l/(4s^ - 1). The SM 
point (Cf M , Cfo 1 ) is marked by the red dot. 



Including the semileptonic decay data at 68% (95%) C.L. we obtain 

RM < 4% (4%), RM < 16% (28%) , (li; 



cutting into the models' parameter space, cf. Eq. (10). Due to the Z-penguin dom- 



inance of the SUSY flavor contributions to Cg and C±o, the constraints Eq. (11) are 



much stronger than those on model-independent scenarios with general Cg and C\q [6]. 
2.4 Constraining (S^lr 

Previous works found that (#23)^ is essentially unconstrained by \AB\ = \AS\ = 1 



decay data 40 . We illustrate the impact of the new constraints from semileptonic 
decays on squark flavor mixing. In Fig. [2]we show the upper limit on | (^23)l_r l> without 
(left-hand plot) and including (right-hand plot) the recent data on semileptonic decays. 
For the latter plot we employ the 68% C.L. bounds that are depicted in dark gray in 
Fig. [TJ For defmiteness, we choose the A t -mi R plane around the SUSY parameter point 
defined in Table [2] One observes that, at least in this part of the SUSY parameter 
space, the bound on (<?>23)l.r is significantly improved by the new data. 

In Fig. [3] we show the upper limits on | (5% 3 ) LR \ in the /i-M 2 and the m^-tan (3 
plane, in both cases including the data from B — > R(*H + l~ . The bounds grow stronger 
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Figure 2: Upper bounds on | (£23)1^ | around the SUSY parameter point defined in 
Table [2] Left-hand side: without semileptonic bounds. Right-hand side: including the 
semileptonic bounds at 68% C.L. 



for decreasing Mi and decreasing mi R . Note that the Higgs mass limits and bounds 
by electroweak precision tests are discarded in Figs. [2] and [3j as we only want to show 
here the parametric dependence of the new semileptonic bounds. 

Note that (523) lr is also constrained by demanding vacuum stability, 



OS 



23)LR 



< 



Ml 



2M 2 U + M 2 



1 ■ 



(12) 



— 2 — 2 

where M u and M l denote an averaged up-squark and slepton mass squared, respec- 
tively ETJ. Requiring only metastability, these bounds are in general weakened but 
only very little so in case of {5% 3 )lr 42 . Depending on the flavor-diagonal MSSM 
parameters, the bounds on (^23)1^ from B — » R(*H + l~ data which we obtain are 



stronger than the vacuum stability bounds. For instance, setting m* 



mj 



1 TeV 



and mi = 300 GeV, vacuum stability requires (#23)1^ ~ 0-3, which is comparable 
to the bounds shown in the left-hand plot of Fig. [2] but weaker than the ones in the 
right-hand plot, which includes the new semileptonic rare decay data. 



3 Predictions 



We present implications for 6-physics in Section 3.1, for RFV models in Section 3.2 



and for rare decays of the top in Section 3.3 



3.1 b- physics 



The constraints Eq. (11) imply upper bounds on CP phases of the same order in 



the respective Wilson coefficients assuming order one phases in the flavor parameter 
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Figure 3: Upper bounds on | (S^) LR | around the SUSY point given in Table [2] varying 
in different directions around the example SUSY parameter point. Both plots include 
the 68% C.L. semileptonic bounds. 



(^2s)lR- Only little is known to date about the CP-violating phases of Cg^: the 
relative phase |argC9C* | is constrained by data on B — > K*l + l~ decays at large 
dilepton masses to be near 7r 43 . The situation will improve in the future with more 



precise data and measurements of CP asymmetries. Specifically the (naive) T-odd CP 
asymmetry (A®) is unsuppressed by strong phases 44 and is directly sensitive to CP 



violation in C±o and hence SUSY flavor. At large dilepton masses the CP asymmetry 

(3). 



2N 



> is promising. The latter is related to the CP asymmetry of the forward- 
backward asymmetry put forward in Ref. 
We predict from Fig. [Tl that the B„ — > 

/I lev 12 



/i 1 /i~ branching ratio B(B S — > fjT ^) oc 
i B |^io | is enhanced (suppressed) with respect to the SM one by at most a factor of 
1.3 (0.5) at 95% C.L. and within the range 

1 x KT 9 < B(B S < 5(6) x KT 9 . (13) 

On the other hand, taking the recent upper limit Eq. ([!]) one obtains at 95% C.L. ^] 

Riofa) < 46% (20%) , (14) 



consistent with Eq. (|ll|). 
constant fs t 

B(B S -)• yU + /i _ ) S M 

f Ba = 256(6) (6) MeV from Ref. |47], with B(B, 



In Eqs. (13) and (14) we used for the 5 s -meson's decay 
231(15)(4) MeV from Ref. [46], corresponding to a SM branching ratio 
(3.1 ± 0.6) x 10~ 9 . The numbers in parentheses correspond to 

/ u + /i")sM = (3.8 ± 0.4) x 10~ 9 . In 



the lattice results the first error is statistical and the second one from systematics. 



2 It has been pointed out recently [45] that the finite lifetime difference in the B s -system causes 
B(B S — > when extracted from an untagged measurement as in Eq. (lip to differ from its 



corresponding value in the unmixed case. Including the effects from mixing, the limits in Eq. (14) 
would get stronger by a few percent. 
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One observes that the purely leptonic and the semileptonic decays give about com- 
parable constraints. Since the former has a larger sensitivity to scalar/pseudoscalar 
operators, which would in SUSY kick in for large values of tan /3 and a not too heavy 
Higgs sector and which we neglect, the combined analysis of both modes is most im- 
portant as they probe complementary NP. 

Furthermore, the forward-backward asymmetry in B — > K*l + l~ decays exhibits a 
well-known zero in the SM, roughly determined by — 2m fe m s Re(CV(Ai&)/C , 9(/i;,)). In 
the SUSY model both C-j and Cg are near their respective SM values, and so is the 
location of the zero. This is consistent with data, cf. Eq. (|2|. 

3.2 Implications for RFV models 
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Figure 4: The required value of (S^lr in the RFV model, see Eq. (15), for m 
GeV and (5* 3 ) LR = 0. 



1000 



As shown in Section 2.4 the current bounds on (S^lr can reach a level of ~ 0.1. 
This implies constraints on SUSY flavor models which have rather largish values of 



($23) lr m this ballpark. One such model 18 19] (see references therein for earlier 
works) is based on radiative flavor violation, where the small fermion masses and CKM 



off-diagonal elements originate from quantum loops 48 . The CKM matrix is assumed 



to be the unit matrix at tree level, and the off-diagonal elements are induced by quan- 
tum corrections involving off-diagonal trilinear SUSY breaking couplings. Specifically, 



the quark mixing between the second and the third generation is then given as 18,19 



V, 



cb 



2a s 
3irmg 



LR 



C (x,x) 



(AS 



23JLR 



m t 



Co(x,y) 



(15) 
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where x = rn~/rin£ and y 



mt~ /ml and 



C (x,y) 



[1 - y)x\og(x) + (x- l)ylog(y) 
(x-l)(y- \){x-y) 



(16) 



The dimensionless loop function Co satisfies Cq(1, 1) = —1/2 and Cq(1,0) = — 1. 
Above, we allowed for an admixture of flavor generation from squark mixing in the 
down-sector (^23)^ = (^23) LJ?/ m l as we ^" 
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Figure 5: Excluded regions of RFV parameter space in the rrig — m- tR plane around 
the SUSY point given in Table [2] with M 2 = 800 GeV. The left-hand plot displays 
all bounds including the semileptonic bound at 95% CL. (black line), the limit from 
m t ~ > 100 GeV (red line), a hypothetical limit m t - > 400 GeV (orange line), (green 
line) and b — > 57 (blue line) for rrig = 1000 GeV and (S23) LR = 0- m the right-hand 
plot the 95% C.L. semileptonic bounds are shown for different values of m g and in 
addition with finite (#23) M = ±1 x 10 -3 



From Eq. (15) one can determine the flavor mixings required for a realistic CKM 
(40.6 ± 1.3) x 10- 3 f 



element V cb = (40.6 ± 1.3) x 10~ 3 [26] . Since C < 0, it follows that (5% 3 ) LR > is 
required, and that an additional contribution from (S^lr > (< 0) demands a larger 
(smaller) value for (S% 3 )lr. The value of (c^) lr required to generate V c b is shown in 



Fig. § for m g = 1000 GeV and (5^) LR = 0. We evaluate Eq. (|15j at fM = 120 GeV. 



Since the requisite (5% 3 ) lr is sufficiently below one, at least for the region displayed and 
for reasonable gluino masses, the use of the mass insertion approximation is justified, 
For the gluino-down-squark contribution to the Wilson coefficients we 



sec also 19 



use the mass insertion approximation of the exact results given in 14 . 

We further evaluate the constraints from K° — ^-mixing via up-squark mixing 
with the third generation, induced by (8% 3 ) LR (5i 3 ) lr- The latter factor has a large CP- 
phase since (Si 3 )lr generates V u b through the analogue of Eq. (15) 19 . We use the 
full, non-MI result given in Ref. [49] with erratum 50 and applye^ F ^< 0.6e^ p r] We 



3 We thank Andreas Crivellin for numerical checks of the chargino contribution to K° — if °-mixing. 
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SM 


MSSM pre-LHC 


ATLAS 10fb _i 


ATLAS 100 fb" 1 


t — > cy 


4.6 x 10" 14 


5.2 x 10" 7 


9.4 x 10" 5 


3.0 x 10" 5 


t -)■ eg 


4.6 x 1(T 12 


3.2 x 10" 5 


4.3 x 10" 3 


1.4 x 10" 3 


t cZ 


1 x 10~ 14 


1.8 x 1(T 6 


4.4 x 10" 4 


1.4 x 10" 4 



Table 5: The branching ratios of rare top decays in the SM 54[ 55 and upper limits 
in the MSSM from the pre-LHC era 



sensitivity (5a observation) with 10 fb 



56 
i 



The last two columns denote the ATLAS 
and 100 fb -1 , 



respectively |57| . 



use the NLO-RG factor 77 ~ 0.8 p31 52 for the leading AF = 2 operator [s7 M (l — 7s)rf] 



2 



53 



with bag parameter B% h (2 GeV) = 0.52 

The interplay of the various constraints is illustrated in Fig. [5j We show the ex- 
clusion regions in the mi R -rriq plane with SUSY parameters chosen around the SUSY 
point defined in Table [2] with M 2 = 800 GeV. Note that in the chargino-loop b — > sl + l~ 
amplitude some CKM elements of the Wilson coefficients, see, e.g., Eq. ^ in the MI 
approximation, cancel against those in Eq. (|3| to only diagonal elements V t bV* s . For 
those we use the physical CKM matrix elements instead of the bare ones [19] . We find 
that the constraints from ex are weaker than the semileptonic bounds for not too small 
values of M 2 and for stops sufficiently split from the other squarks, around the point Ta- 
ble^ for M 2 > 500 GeV and m^ R < 800 GeV. In other words, the Glashow-Iliopoulos- 
Maiani (GIM) suppression can efficiently make evade the semileptonic bounds. We 
learn that the new B — > K^l + l~ data allows to exclude an additional part of the pa- 
rameter space. The excluded region is larger for lighter gluinos and for contributions 
from {5i 3 ) LR > 0. 

3.3 Implications for rare top decays 

Since the same flavor- changing squark mass parameters enter B-meson- and top quark 
decay amplitudes, the improved constraints on (^3) LR obtained in Section [2] can lead 
to sharper upper bounds on the branching ratio of rare top decays in the MSSM. 
Specifically, we consider t — > cV, where V is either a photon, gluon or a Z-boson. In 



the SM, the branching ratios are negligibly small, related to a GIM suppression 54 55 
Estimates together with upper limits from a pre-LHC study for the MSSM 56 and 
the requisite branching ratios for 5cr observations at center-of-mass energy -y/s = 14 
TeV with 10 fb -1 and 100 fb" 1 of collision data [57] are compiled in Table [5J Thus, an 
observation of these decays at the LHC is already excluded for the case of the MSSM. 
In order to further strengthen this point, we nevertheless find it useful to provide 
updated bounds on the achievable branching ratios in the MSSM. 

The leading contributions arise from squark-gluino loops with non-vanishing (S^) lr 
and/or [S^ll- It turns out that the largest effect in the rate into all three final states 



is due to diagrams involving (5% 3 )lr, cf. 56 . We calculate the branching ratios for 



each parameter point in the region defined by Table |4j which passes the constraints i) 
to v) as well as the constraints on Cg and C\q from B — > R(*H + l~. The gluino mass 
is taken to be rrig = 700 GeV. For heavier gluinos the upper bounds on the branching 
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ratios decrease further. We employ the formulae for the t — > cV decay widths which 



are found in 58 . The branching ratios are obtained by normalizing F(t — >■ cV) to the 



SM value for the dominant decay mode of the top quark, T(t — > bW) = 1.29 GeV 26 
Note that, as in Section [2j we do not rely on the MI approximation. 

As a result, we find that the maximal branching ratios compatible with the 68% 
C.L. constraints from the rare -B-decays are 

®{t cj) < 2.1 x HT 8 , B(t eg) < 7.2 x 1(T 7 , B(t cZ) < 1.0 x 1(T 7 . (17) 

The constraints improve on the ones from a previous pre-LHC study shown in Table 
[5] by a factor of 24, 44 and 18 for decays to 7, g and Z plus charm, respectively. The 
improvement is, depending on the parameter point, due to the AB = 1 data but also 
from the improved squark and gluino mass bounds from the LHC Without the Higgs 



mass constraint iv) enforced using FeynHiggs v.2.9.0-beta, the limits in Eq. (17) would 
be lifted by about 0(5 - 8). 



4 Summary 

Recent progress in both the theoretical description and the data regarding B — > 
K^l + l~ decay distributions allows to put new constraints on squark flavor viola- 
tion. While the current sensitivity in C\q is about a factor of two away from MFV 
physics, beyond-CKM flavor mixing is strongly constrained, see Fig. [1} This concerns 
prominently left-right flavor mixing between the second and the third generation in the 
up-sector encoded in the parameter (S^lr- With its sensitivity in other rare processes 
including B s — B s mixing and b — > S7 modes being very weak, (#23) lr has previously 
been bounded very loosely only. 

We obtain - depending on the flavor-diagonal SUSY parameters - constraints as 
low as (#23) lr ^0.1. This excludes solutions to the flavor problem with RFV flavor 
models based on flavor generation in the up-sector and sub-TeV spectra. Models with 



horizontal flavor symmetries generically predict (5% 3 )lr ~ V c b(m t /mq), cf. e.g., 59 
and are about an order of magnitude below the current limits. Note that in the future 
(#23) lr, could be probed in the presence of an appreciable first-third generation mixing 



in K — 7- ixvv decays |60| at the NA62 experiment 61 



The new flavor constraints lead to predictions for b- and top physics: 

• The B s — > branching ratio is, assuming no scalar /pseudoscalar contribu- 



tions, bounded from below, at about ~ 1 x 10 9 , see Eq. (13). The constraint 
on the short-distance coupling C±o from the experimental upper limit Eq. Q is 
consistent with the one obtained from semileptonic decay data. Both searches 
should be pursued further as they probe complementary NP. 

The forward-backward asymmetry in B — > K*l + l~ decays has a zero at low 
dilepton mass. The position is near its SM value, consistent with the preliminary 
determination by the LHCb collaboration, see Eq. (|2]). The forward-backward 
asymmetry of the recently observed B s — > \T and A b — > Afi + [i~ decays 



shares the same features 62 
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• Top FCNCs t — > cj,t — > cg,t — > cZ have branching ratios below < few x 
10~ 8 , 10~ 6 and 10 -7 , respectively, and are too small to be observed at the LHC 
with foreseeable luminosities. 

The | AS | = | AS | = 1 constraints will improve in the near future with the large 
b — > sZ + /~-induced event samples expected from the LHC experiments. Progress arises 
from better statistics in combination with the availability of additional observables. We 
highlight here those sensitive to CP violation in Cio - induced by a complex-valued 

• The (naive) T-odd CP asymmetry (Af) |44j can be up to order ten percent at 
small dilepton masses [13] . 

• The CP asymmetry (a^p) does not require flavor tagging. At large dilepton 



masses it can reach a few percent 43 . 



The flavor constraints reported here get stronger for lighter stops, an ingredient of 



TeV-scale model building, recently, e.g., 63-66 , with possibilities to be seen directly 



at the LHC. If realized in nature, chances are that such models or others based on the 
generic MSSM with light stops will show up one way or the other, or both. 



Note added: During the publishing process a related and complementary work 67 
on MFV SUSY constraints at large tan f3 from rare S-decay data appeared. 
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A Loop functions 

We derived the dependence of CV^io on (5 23 )lr in the MI approximation, cf. Eq. (J7|, 
from the exact expressions in 14 . We find the following loop functions: 



F(x 1 ,x 2 ,x 



t R j 



* 2 fl(Xi/x iR ) - fi(Xi) 



i= l2 Xi/Xl R -Xi 



(A.l) 



F Z - p (x 1} X 2} Xt R ) = X av ^ V jl V i2 {^ii^iiV ■' '■'./ 

»J=1,2 



C 2 (xx ,Xi,Xj) -c 2 (l,x i ,x j ) 



+25, 



C 2 (Xj,l,X iR ) -C 2 (Xj, 1,1) 



% - 1 



F 7 " p (xi,x 2 ,x 4fl ) 



where 

27 j 1 j X{ 



9' 



X 



tR 



t=l,2 



=4z <M £V il V£|V$ 1 | 



Xj/x iR f 7 {Xi/x iR ) - Xif 7 {Xj) 

d 2 (x j, Xj , x^ 5 xpj ) d 2 {xi , , 1 , Xj~, 1 ) 



(A.2) 
(A.3) 

,(A.4) 



m 



m 



tR 



7< 



(A.5) 



and x a „ = |(5 + xz). The functions /j, Cj, d 2 are defined in 14 . The unitary matrices 



U and V which appear above diagonalize the chargino mass matrix: 

U*M^±V ] = diag (m^,m^ 2 ) , 

where 



M2 \/2mw sin /? 
\f2mw cos /3 /i 



(A.6) 



(A.7) 



The up-squark mass matrix is parametrized as 



M 



2 = fm) LL (Mi) LR 
5 V(M|)L (M?) RR 



(A.8) 



with (M|) 



'0 

( M *)lk = (00(A« 3 ) M ) , (M s 2 ) 



v (A&) Lfl/ 



RR 



'm\ 
m\ 
m? 



(A.9) 
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with the diagonal elements m~, (A^)^ = mt (A t — //cot 0) and to| , defined as 

2 

m L = m L,soft + m t + 3 S w m l cos 2/3, (A.10) 

where m?^ goft denotes the pure soft term contribution. The stop masses are obtained 
after diagonalization of Mj: 




< + m l^v H - m l) 2 + 4 ( A &)L + 4 ( A 3 3 )L J • (a- 11 ) 



The other squark masses are then given as m„-. 



Our findings for the Wilson coefficients CV^io agree with Ref. 16 if the following 



modifications are made to the results of Ref. 16 



1. ) Eqs. (25) and (26), (7-penguin with chargino loop, contributions to Cg): 1/3P 42 

has to be replaced by — 1/3P 42- 

2. ) Eq. (26), (7-penguin with chargino loop, contributions to C 7 ): The complex 

conjugation of has to be removed. 

3. ) Eq. (31), (7-penguin with gluino loop, contributions to C 7 ): The (5$ 3 )ll contri- 

bution has to be multiplied by 2 and (S$ 3 )rl has to be replaced by — (5$ 3 )lr- 

4. ) Eq. (34), (7-penguin with chargino loop, case of light stop, contributions to Cg): 

The expression has to be multiplied by a factor 1/6. 

We explicitly checked the above corrections by direct calculations using the conventions 
of 14 . Note that the relative minus sign between 16 and 14 in front of (5$ 3 )lr m the 



7-penguin contribution to C 7 , see item 3., can be absorbed by a rephasing of the gluino 
mass and the A terms, or equivalently, by rephasing the right-chiral down squarks in 



the super-CKM basis. The corrections 1., 2. and 4. are in agreement with Ref. 68 . 
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